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Introduction 

Chronopotentiowitry  <1)  has  not  ba«n  a  widaly  used  alactroanalytlcal 
tachnlqua  dua  to  tha  savara  distortion  of  the  predicted  shape  of  the  response 
due  to  capacitativa  charging  of  tha  electrode.  Vhile  tha  aagnitude  of  the 
total  currant  applied  to  tha  electrode  can  be  controlled,  tha  fraction  of  the 
total  currant  driving  the  faradalc  reaction  is  tiae  dependant.  As  tha 
electrode  potential  changes  as  tha  logaritha  of  the  ratio  of  the  surface 
activities  of  tha  redox  couple,  so  Bust  the  the  aagnitude  of  the  capacitativa 
charging  current .  .This  current  depends  on  the  rate  of  potential  change 


according  to  C  (dE/ 


x  auaa 

E/At) . 


where  C  is  the  double  layer  differential 

U 


capacitance.  When  the  reactant  surface  concentration  approaches  zero  (the 
transition  tiae),  the  rate  of  change  of  potential  is  high,  and  the  distortion 
of  the  predicted  shape  tends  to  Increase.  The  tiae  to  the  transition  is 
important  in  electrocheaical  analysis  and  we  find  that  it  can  therefore  not  be 
aeasured  easily  with  a  high  degree  of  accuracy.  Elaborate  instruaentatlon  has 
been  designed  to  offset  the  capacitative  current;  these  Include  feedback 
compensation  circuits  (2)  and  the  use  of  blank  cells  (3). 

v^t  is  well  known  that  one  of  the  aajor  advantages  of  aicroelectrodes  (4) 
is  the  reduced  aagnitude  of  the  capacitative  component  of  the  current  (which 
varies  as  the  area  of  the  electrode)  with  respect  to  the  faradaic  component 
(which  usually  varies  as  a  one  dimensional  characteristic  length  of  the 
electrode,  such  as  the  radius  of  a  disk).  These  devices  should  therefore  be 

ideally  suited  to  galvanostatic  measurements ^  Aoki  et  al  (5)  have  presented 

\ 

an  approximate  expression  for  the  response  at  disk  electrodes  which  is  based 
only  on  a  constant  surface  concentration  model  whereas  a  constant  flux  model 

is  probably  more  applicable  (6)  to  microelectrodes;  the  exact  solution  lying 

! 

somewhere  in  between  these  two  limiting  models.  In  any  case,  the  results  for 
both  cases  are  similar. 


OK 
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In  our  work  (6) ,  wo  reported  the  exact  integrals  of  the  time  dependent 
diffusion  equations  dascribing  aass  transport  to  disk  electrodes  of  finite 
size.  The  result  follows  froa  a  general  analysis  based  on  the  properties  of 
discontinuous  integrals  of  Bessel  functions.  The  aethod  has  been  used  to 
solve  a  variety  of  probleas  in  heat  conduction  (7)  as  well  as  electrochemistry 
(4,8-11).  In  this  paper  we  describe  the  fora  of  the  chronopotentlograa  for 
constant  laposed  flux  conditions,  and  the  analogous  response  to  a  linearly 
swept  current  experiment. 

Kxperlaental 

The  instrumentation  consisted  of  a  low-noise  galvanostat  built  in  these 
laboratories,  and  a  Hl-Tak  Instruments  nodal  PPR1  waveform  generator.  The 
working  electrode  was  tied  to  ground  through  a  resistor,  and  the  potential  of 
the  working  electrode  was  measured  with  respect  to  the  reference  with  a  high 
impedance  differential  aaplifier.  The  output  was  measured  with  an  x-y 
plotter,  or  with  a  Lecroy  model  9400  digitizing  oscilloscope.  The 
electrochemical  cell  consisted  of  a  large  platinum  wire  secondary  electrode 
palladium  sheet  reference  electrode,  and  a  disk  microelectrode.  The  cell  was 
purged  with  nitrogen  before  measurements.  The  entire  cell  assembly  was 
rigidly  mounted  Inside  an  aluminum  Faraday  cage  (to  reduce  capacitatively 
coupled  noise)  which  contained  triaxial  bulkhead  connectors  for  feeding 
through  the  instrumentation  cables. 

The  disk  microelectrodes  were  prepared  by  heat  sealing  the  appropriate 
fine  wires  into  glass  tubes.  The  disks  were  polished  with  fine  alumina  (0.3 
and  0.05jim  particle  diameters,  Baikalox  International)  rinsed  thoroughly,  and 
sonicated  briefly  before  use.  Several  disk  sizes  from  1.0  to  2S.0pm  were  used 
in  the  experiments. 

Aqueous  solutions  were  prepared  using  analytical  grade  reagents  in  18Mfl 
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Nanopur*  v«C*r.  Potassium  ferrocyanld*  (Mallinckrodt) ,  potassium  ferrlcyanide , 
potassium  ehlorlda,  and  potassium  nitrate  (J.  T.  Baker)  were  used  as  received. 
The  ferri-  and  ferrocyanide  were  present  in  equimolar  concentrations  in  0.1M 
potassium  nitrate.  Rutheniun(IlI)  hexamine  chloride  was  prepared  and 
purified  froa  the  chloride  by  standard  procedures. 

In  the  linear  sweep  aaperoaetry  experiments,  the  sweep  rate  was  varied  up 
to  50pA  s'1;  the  starting  value  was  varied  froa  -lOnA  to  OnA,  and  the  ending 
value  was  varied  froa  -lOnA  to  +10nA. 

Discussion 

For  a  simple  electrochemical  experiment  involving  a  single  reactant,  w* 
have  solved  the  time  dependent  diffusion  equation 


3C  p  a2c  ,  D  8C  t  p  d2C 
a t  8x2  r  ar  aza 


[1] 


in  circular  cylindrical  coordinates  (6).  Here  C  is  the  concentration  of  the 
reactant,  and  r  is  the  radial  distance  coordinate  measured  froa  the  center  of 
the  disk  electrode  which  is  imbedded  in  the  insulating  plane  at  z  ■»  0.  The 
general  initial  condition  is 

r  >  0,  z  >  0,  t  -  0,  C  -  C*  [2] 

W 

where  C  is  the  bulk  concentration.  For  the  chronopo tent iome try  case,  we  take 
a  constant  uniform  flux  -Q  (mols  cm*2s*1)  over  the  surface  at  all  times  t  >  0. 


The  boundary  conditions  at  the  surface  of  the  electrode  are  then 
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0  i  r  <  .,  .  -  0.  t  >  0  .  d[!§]  -  -  q. 

r  >  -  0,  t  >  0.  »[s]  -  »• 

from  which  we  derived  the  integral  for  Equation  [1] 

C  "  C"'  D5  f  J„<«) J^oajerf <0X '  2otX/ *)^[ 

o 

tfe  found  that  the  average  surface  concentration  is 


C*r  *  C"‘  ^  J  [v«>]a*rf<D1/2«tl/a)^! 

With  the  substitutions 


[3] 


[A] 


(51 


l2  -  Dt 
fi  -  al 

ve  can  write  [5]  in  the  aore  useful  dlnensionless  form 


i2 


(6] 

(7) 


(8] 


The  function  has  been  tabulated  (6)  as  a  function  of  the  dimensionless 
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parameter  (Dt/e*) ,  but  la  readily  evaluated  by  appropriate  numerical 
Integration  techniques  (excellent  accuracy  is  obtained  using  a  modified 
Burlirsch-Stoar  method  (12)). 

At  long  times  and  sufficiently  small  values  of  the  flux,  a  transition 
time  is  not  observed,  and  we  have  the  steady  state  condition  (6) 


If  Q  Is  sufficiently  large  we  observe  a  sharp  transition  as  the  surface 
concentration  of  the  reactant  approaches  zero,  and  we  have  the  result 


from  which  we  can  calculate  the  transition  time. 

The  requirement  on  2Qa/DC*  for  steady  state  conditions  has  been  shown  (6) 

to  be 

^  S  ~  •  2.35619  [11] 

4 

DC 


For  the  case  of  simple  Butler-Volmer  kinetics  and  for  equal  concentrations  of 
both  the  oxidized  and  reduced  forms  of  the  redox  couple 
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Q 


•xp 


(l-a)nF 


Rt 


I 


[12] 


where  we  use  the  average  concentration,  Equation  [8].  Ve  obtain  an  expression 
[13]  from  which  the  polarization  plots  aay  be  determined: 


*9  +  £&  • 

o  DC 


Dt 

l 


jexp 


RT 


+  exp 


(l-q)qF 


RT 


'  ' 

t 

■  ► 

4  4 


-  exp 


-arjF 

RT 


-  exp 


(1“«)«»F 


RT 


U3] 


The  applicability  of  the  use  of  the  average  surface  concentration  for 
aicroelectrode  experiment*  has  been  discussed  in  detail  (6).  This  expression 
shows  that  the  transients  are  a  function  of  a,  FQ/i  and  2Qa/Dc"\  It  is  clear 

P 

that  the  parameter  FQ/i  determines  the  reversibility  of  the  reaction,  while 

O 

2Qa/DC*  controls  the  ratio  of  the  imposed  reaction  rate  Q  to  the  mass  transfer 
rate.  Figures  1*3  show  the  shapes  of  the  chronopotentlograms  as  a  function  of 

2  co 

the  dimensionless  time  variable  Dt/a  and  the  parameters  PQ/i  and  2Qa/DC  . 

O 

The  roots  to  Equation  [13]  were  obtained  by  Brent's  method  (12).  Ve  point  out 
once  more  that  for  low  values  of  the  mass  transport  parameter,  Equation  [11], 
there  is  no  transition  (Figures  1*3).  This  is  due  to  the  fast  approach  to 
steady  state  diffusion  at  small  disk  electrodes.  It  is  therefore  probable 
that  constant  current  experiments  are  not  particularly  useful  for  analytical 
applications,  since  it  will  be  difficult  to  choose  an  appropriate  flux  for 
observation  of  a  transition  time.  A  simple  linearly  swept  current  of  the  form 
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Q(t)  -  7t  (14] 

however,  will  always  give  rise  to  a  transition  sines  as  Q  increases,  2Qa/DC* 
must  eventually  exceed  3*/4.  By  a  siailar  analysis  (6)  we  have  found  the 
relation,  for  the  Butler  Volaer  kinetic  nodal  (Equation  [12]) 


Fyt  4yta 
io  DC* 


-  **P 


isnl 

RT 


-  «P 


(l-o)f>F 


RT 


[15] 


where  is  given  by 


y  erf (y)  dy 


d£ 


[16] 


which  also  has  been  tabulated  (6)  and  can  be  evaluated  similarly  to 

The  roots  to  Equation  [15]  were  again  found  by  Brent's  method  (12).  In  this 

case,  we  see  that  yt  is  not  independent  of  the  dimensionless  time  variable 
2 

Dt/a  .  We  have,  however, 


4yta  4ya3 
DC*  "  D2C* 


[17] 


and 


He  have  therefore  evelueted  the  overpotential  i|  n  i  function  of  the  verieble 
Dt/a2  and  of  the  parameters  FyaVi  and  47a3/  DlC*.  The  results  are  presented 

O 

in  Figures  4*6,  which  are  for  faat,  intermediate ,  and  slow  rates  of  reaction, 
respectively.  Unlike  the  constant  current  results  shown  in  Figures  1*3,  a 
transition  is  always  observed. 

Figure  7  shows  the  sharp  transition  that  is  observed  in  a  typical 

experiment.  The  response  here  is  for  a  linear  sweep  amperoaetry  experiment  at. 

a  platinum  disk  microelectrode .  The  current  was  initiated  at  a  large  negative  > 

value,  and  swept  towards  positive  values.  This  result  is  in  contrast  to  the 

poorly  resolved  transitions  obtained  with  electrodes  of  conventional  size 

where  there  are  large  capacitative  effects. 

Figures  8-9  are  responses  for  linearly  swept  current  experiments  in  which 

the  initial  current  was  zero;  these  experiments  can  be  therefore  be  analyzed 

directly  by  Equation  [IS].  The  size  of  the  disk  electrode  was  determined  by 

measuring  the  diameter  of  the  wire  before  mounting  in  its  holder,  by  optical 

microscopy,  and  by  measuring  the  limiting  currents  in  standard  redox 

solutions.  These  methods  gave  results  that  agreed  to  within  4%.  The 

2 

responses  were  made  dimensionless  and  plotted  vs.  Dt/a  ,  and  appear  in  Figures 
8  and  9  along  with  the  closest  fit  theoretical  plot.  The  data  was  converted 
by  noting  that 

Q(t)  -  yt  -  l(t)/F  -  i(t)/AF  [19] 

where  A  is  the  surface  area  of  the  disk,  and  I(t)  the  imposed  current  density 


at  time  t.  If  w  i i  the  current  sweep  rate  (A/s),  then  Equations  [17]  and  [18] 
can  written  in  the  fora 


4yt 

DC* 


’  ' 

Dt 

4  va3 

r  ' 

Dt 

4va 

*  ' 

Dt 

2 

*  J 

d2c*af 

[«2j 

wD2c"F 

2 

* 

[20] 


and 


Fjt 

1 


Ft/aa  Dt 

DAFi  *  a 

O  A 


[21] 


Ve  have  assumed  a  value  for  the  diffuaion  coefficient  of  6.7  x  10*#cm*s'1 
for  both  redox  species,  and  a  value  of  the  exchange  current  density  i  of  9.0 

O 

wA/cm2  for  the  equiaolar  laK  solution. The  paraaeters  were  then  used  to 
calculate  the  theoretical  responses  shown  as  the  solid  lines  in  Figures  8  and 
9.  For  a  fixed  sweep  rate,  the  increase  in  the  current  density  with  decrease 
of  electrode  size  leads  to  an  increasingly  rapid  transition.  For  the  smallest 
electrode  used  (lpm  diameter)  the  transition  (when  using  the  saae  current 
sweep  rate  v  as  for  the  other  electrodes)  was  very  rapid  as  expected,  but 
analysis  was  not  possible  due  to  the  finite  response  times  of  the  potential 
measuring  equipment. 

The  potential- time  response  for  platinum  microelectrodes  under  other 
imposed  flux  conditions  is  shown  in  Figure  10.  In  these  experiments,  the 
current  sweep  was  initiated  at  negative  values,  and  the  responses  are 
therefore  not  predicted  directly  by  the  above  analysis.  The  qualitatively 
expected  behavior,  however,  is  clearly  evident  in  the  Figures.  At  the  larger 
(25pm)  electrodes  steady  state  diffusion  is  not  developed  rapidly.  For  a 
given  value  of  the  initial  applied  flux  and  sweep  rate  one  therefore  expects 


the  transition  to  positive  potentials  to  be  delayed  by  the  non* steady  state 
diffusion  of  the  two  redox  species.  The  oxidation  of  the  excess  ferrocyanide 
(generated  by  the  initial  negative  flux)  clearly  inhibits  the  tiae  to 
transition  to  positive  potentials.  Figure  10a;  in  fact  we  see  that  the 
potential  is  substantially  negative  when  the  laposed  flux  la  zero.  The  effect 
is  considerably  reduced  at  the  lOpa  electrode,  and  is  not  discernible  at  the 
lim  electrodes,  Figures  10b  and  c,  at  which  the  aass  transport  rates  are 
significantly  enhanced. 

The  use  of  disk  aicroelectrodes  of  suitable  radii  in  galvanostatic 
experiaenta  therefore  appears  to  yield  vary  predictable  and  diagnostic 
analytical  results.  The  sharp  transitions  observed  in  linear  sweep 
aaperoaetry  experiments  assures  accurate  aeasureaent  of  the  transition  times 
and  potentials.  The  use  of  aora  complex  current  sweep  programs  has  other 
advantages  in  the  study  of  coupled  chemical  reactions,  and  will  be  the  subject 
of  a  further  report  shortly. 
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Glossary  of  Symbols  Ussd 


A 

C 

C* 


41 

D 

K 

P 

Kt) 

ICt) 

lo 

J.*Ji 

Q 

Q(t) 

R 


r 

t 

T 


y 

z 


a 

a 


7 

i 


Radius  of  disk,  cn 

Elactroda  araa,  cn*2 

Concentration,  sols  cn*9 

Bulk  concentration,  sols  cn  9 

Avarags  concentration,  sols  cn'9 

Double  layer  capacitance,  C/cna 

Diffusion  coefficient,  cn2s  1 

Electrode  potential,  V 

Faraday  constant,  96485  C  equivalent*1 

Current  at  tine  t,  A 

Current  density  at  tine  t,  A  cn*2 

Exchange  current  density,  A  cn*2 

Bessel  functions 
*2  "1 

Flux,  sols  cn  s 

•  2  *1 

Flux  at  tine  t,  nols  cn  s 
Gas  constant,  6.314  J  molsHc  1 
Radial  coordinate ,  cn 
Tine,  s 

Temperature,  K 

Dunny  Integration  variable. 

Coordinate  normal  to  plane  of  disk,  cn 
Transfer  coefficient  (when  in  exponent) 
Continuous  dumny  integration  variable. 
Parameter  ai 

-  j  -2 

Flux  sweep  rate,  nols  cm  s 
(Dt)1/l 
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ij  Ovarpotantlal ,  V 

v  Currant  avaap  rata,  A  a 


IS 

Copy  available  to  DTIC  dee#  not 

Figure  legends  p®nnit  ^  J«<nl>Ie  reproduction 

1.  Chronopo tent iogr  sms  as  predicted  by  Equation  13  for  totally  revere lble 
electrode  kinetica  (1  -  •>,  FQ/1  -  0)  as  a  function  of  the  aass  transport 

O  O 

paraaeter  2Qa/DC  and  the  diaensionless  tine  variable  Dt/a  . 

2.  As  Figure  1  except  for  intermediate  values  of  the  kinetic  paraaeter 
(PQ/i-  1.0,  quae irreversible  reaction). 

3.  As  Figure  1  except  for  large  values  of  the  kinetic  paraaeter  (FQ/i-  100.0 
-irreversible  reaction). 

i 

4.  Linear  sweep  aaperoaetry  response  as  as  predicted  by  Equation  IS  for 
reversible  electrode  kinetics  (i  -  •,  Fya^/Ll  -  0)  as  a  function  of  the 

O  O 

3  2  oo 

aass  transport  paraaeter  4ya  /D  C  and  the  diaensionless  tlae  variable 
Dt/a2 . 

5.  As  Figure  3  except  for  interaedlate  values  of  the  kinetic  parameter 
(Fya2/Di  -  1.0,  quasi  reversible  kinetics). 

O 

6.  As  Figure  3  except  for  large  values  of  the  kinetic  parameter 
(Fya2/Di  -  100,0,  irreversible  kinetics). 

O 

7.  Example  of  linear  sweep  amperometry.  9.0ms  diameter  gold  disk 
microelectrode,  solution  3.5mK  Ru(NH)3+  in  0.1M  KC1.  Sweep  rate  was 
625pA  a  l,  initiated  from  -l.OnA. 
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8.  Linear  sweep  amperometry  of  lmM  each  Fa(ZI) (CN)**/Fe(III) (CN)'3  in  0.1M 

P  P 

KNO  at  a  25.0pm  dlaaatar  platlnua  microelectrode .  Clrclaa  are 

9 

experimental  data,  and  the  line  ia  the  calculated  response  using  the 
Indicated  parameters.  Current  sweep  rate  v  was  50pA  a"*  Sweep  initiated 
at  O.OOpA. 

9.  Same  as  Figure  8,  except  at  10.0pm  (open  circles)  and  1.0pm  diameter 
(marked  circles)  electrodes. 

10.  Same  as  Figure  8  and  9  except  sweep  Initiated  at  -lOnA. 

-  25pm  diameter  disk 

■  10pm  dlaawter  disk 

.....  l (M  diameter  disk 


OVERPOTENTIAL  /  /  V 


POTENTIAL  /  V  vs.  P* 
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